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Q • Abstract. Statistics of molecular random walks in a fluid is considered with the 

help of Bogolyubov equation for generating functional of distribution functions. An 
invariance group of this equation is found. It results in many exact relations between 
path probability distribution of a test particle and its correlations with the fluid. As 
the consequence, significant restrictions on possible shape of the path distribution do 
' arise. In particular, the hypothetical Gaussian form of long-range asymptotic proves 

to be forbidden, even (and first of all) under the Boltzmann-Grad limit. An allowed 
diffusive asymptotic possesses power-law long tail (cut off by free flight length). 
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1. Introduction 



Random wandering of particles of the matter is mechanism of diffusion and many 



other transport processes as well as various noises and fluctuations. What can be its 



statistics? Strangely, this important question never was addressed to rigorous statistical 

m 
oo 

CO 

lO ' the Lorentz gas asymptotic of random walks is Gaussian [1] , all the more it must be 



mechanics. At present, seemingly, the answer is obvious without it: since even in 



the same in a usual fluid. However, one can notice that such conclusion is founded 

OO 

on instinctive identifying dynamical independence of events in "real life" (concrete 
realization of many-particle system) and statistical independence of events in "theory" 
(statistical ensemble of systems). But already Krylov [2] thoroughly explained that the 
first does not imply the second l. 
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1 In the probability theory by definition A and B are independent if V{A A B) = V{A)V{B) [3]. 
Therefore statement "A and B are independent since have nothing to do with each other" 0] is 
wrong: such A and B may both depend on common conditions so that V(A A B) ^ V(A)V(B) ! 
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This can be understood by the example of gas of iV hard balls in a box if considering 
their motion as motion of one ball in 3iV-dimensional billiard [5] which resembles 
the Lorentz gas. Here as well initially one-dimensional (straight-line) motion becomes 
multidimensional. However hardly it is possible to speak about ergodicity before all the 
dimensions and scatterers become apparent, that is before time ~ |iV(JV— 1) / (N/ r) ~ 
Nt ~ VL/txtIvq passed after start of observations, with Q , r ~ (^r^z/wo) -1 , r , vq 
and v = N/Q being volume of the box, typical free path time of gas particles, radius 
of their interaction (diameter of the balls), their characteristic velocity and mean gas 
density, respectively. Even for 1 cm 3 of the air this is time greater than 1000 years! 
Thus, it was true remark [6] that role of ergodicity in physics is strongly exaggerated 
since in physical reality limit iV — > oo precedes limit t/r — ► oo . 

If it is so then trajectory of test particle in gas of sufficiently many particles (formally, 
under iV — > oo ) is non-ergodic: its even very relatively distant fragments (separated by 
time intervals ^> r ) are statistically dependent one on another although independent 
in dynamical sense. This means that any concrete realization of random walk of the 
test particle (corresponding to some concrete phase trajectory of the whole system) 
possesses its own specific kinetic characteristics (diffusivity and mobility, etc.) or, 
better to say, has no certain kinetic characteristics [8]. 

Theoretical tools for investigations of such statistics a long time ago were presented 
by N. Bogolyubov [9] and his followers. That are the BBGKY (Bogolyubov-Born- 
Green-Kirkwood-Yvon) hierarchy of equations or equivalent Bogolyubov's equation 
for generating functional of many-particle distribution functions. But, unfortunately, 
nobody has learned honest use of them without some cutting off the hierarchy and 
justifying such vivisection by imaginary intuitive "independencies". Therefore comple- 
mentary tools are not unnecessary. 
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In the present work we describe an invariance group of the Bogolyubov equation and 
its solutions. For simplicity, it is formulated directly with reference to the problem 
about random walk of test particle in thermodynamically equilibrium fluid. We want 
to demonstrate that exact relations of the group point to essentially non-Gaussian 
statistical features of long-range asymptotic of the random walk. Importantly, in case 
of gas these relations and their consequences, along with the mentioned characteristic 
time scale ~ Nr , are indifferent to gas density, even in the limit of dilute gas (or the 
Boltzmann-Grad limit). 



2. Equations of molecular random walk 

Let box with volume Vt contains N ^> 1 identical atoms plus one more test particle. 
Atoms have mass m, coordinates Vj and momenta Pj (j = 1,2... N) and interact 
with each other via potential & a ( r j — r k) ■ The test particle has mass M , coordinate 
R, momentum P and interacts with atoms via potential $&(rj — R) . The potentials 
are spherically symmetric and short-range with impenetrable point core. Because of 
interactions the test particle is in chaotic motion, therefore let us name it "molecular 
Brownian particle" (BP). 

We are interested in probability distribution of current position of BP, R(£) , un- 
der condition that at initial time moment t — it was placed at certainly known 
position: R(0) = R , while personal positions of atoms all the times are unknown. 
The simplest statistical ensemble what satisfies this requirement is determined by the 
Liouville equation, dD^/dt = [Hn,Dn], for full normalized distribution function 
of the system, D N , and initial condition 

<J(R-Ro) e- HN ' T 
DN(t = ° ) = JdRjdPj 1 ...J N 5(R-R )e~^ (2 ' 1} 
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to it, where is full Hamiltonian of the system (including interactions with the box 

walls) and f, ... = j j ... dr^dp^ . Evidently, such ensemble differs from the Gibbs 

canonic ensemble by initial BP's localization only. The latter does not prevent us to 

introduce particular distribution functions (DF) F n (t) = Q n f 1 ... f N D N (t) and 

then go to the thermodynamical limit (N oc , Q — > oc , v = N/Q = const ) just 

as in [9]. The only, non-principal, difference from [9J is that all n-atom DF under 

consideration include also BP's variables, i.e. in fact are (n+ l)-particle DF, therefore 

their numeration takes beginning at zero, and in respect to BP's variables all they 

are normalized in usual sense of the probability theory. Writing out complete list of 

arguments, F n (t) = F n (t, R, r^, P, p (n) | R ; u) , where = {rx...r n }, pW = 

{pi... Pn } • In particular, F Q (t, R, P| R ; v ) describes BP itself, and f dH f dP F Q = 

1 . What is for the atomic coordinates, in respect to them all DF are "normalized to 

volume", exactly as in In other words, the "weakening of correlations" between 

distant particles holds: F n — > F n _i G m (pk) , if fc-th particle keeps away from others, 

where G m (p) = (27rTm) _3//2 exp (— p 2 /2Tm) is the Maxwell momentum distribution 

of a particle with mass m . The full Liouville equation induces the BBGKY equations 

dFn r „ „ , d f , „ A d 

dt 



[H n ,F n ] + u— $ 6 '(R-r n+1 )F n+1 + j— / ^(r,- - r n+1 ) F n 

°" Jn+l - =1 °~Pj Jn+1 



(2.2) 

n = 0, 1, . . . ) with initial conditions 



Fn (t = 0) = 5(R - R ) Fi e<?) (r (n) |R; v) G M {P) U]=i G m (Pj) , (2-3) 

where H n is Hamiltonian of subsystem "n atoms plus BP " , $^ b (r) = V$ a ,fe(r) , and 
Fn £q \r^ |R; v) are usual thermodynamically equilibrium DF for n atom in presence 
of BP occupying point R. In principle, that will do for finding F (t, R, P| Ro ; v ) 
and thus probability distribution of BP's path, AR(t) = R(i) — R , without any 
additional assumptions. 



Following Bogolyubov, let us combine all our DF into generating functional (GF) 



n=l " 1 Jn j=l 

and equations (I2.2p into corresponding "generating equation" for it: 

d Z_ L d JL r ( i —\ r 

dt M ' <9R V^' 6^) ' 

where operator £ is composed by usual and variational derivatives, 



(2.4) 



(2.5) 



if N Pi 9 5 

i m or i oip[xi) 



(2.6) 



$ 6 (R-ri), 
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+ 



+ ^// [l + ^i)][l + ^ 2 ) 



5^(xi) 
$ ffl (ri - r 2 ) , 



with Xj = {rj,pj} . This is direct analogue of equation (7.9) from [9j. To make it 
well visible, notice that ip(x) = u(x)/v , where u(x) is functional argument used in [9], 
and factor [l+^sci)] [l + ^(x 2 )] can be replaced by [if)(xx) + ip(x 2 ) + ip{xi)i(}(x2) ] due 
to identity J ± L [$ a ( r i — r 2) , ••• ] = . Initial condition to equation (12.51) is obvious: 

^{0, R, P, | R ; u} = 5(R - R ) G M (P) F {eq) {<P\ R; u} 

oo n 

^){0|R;,} = l+£^ 



(2.7) 



n=l 



where we introduced new functional argument 0(r) = J ^(r, p) G m (p) dp and be- 
sides generating functional T^ eq) of equilibrium DF. It is easy to guess that expression 
Gm(P) ^' e9 ' should bring stationary solution of (12.51) : 



P/M) • «9/<9R + £ 1 G M (P)^ eq) {(P\R;u} = 



From here equation 



d $/(r - R) 



<9r 



+ 



T 



SJT(eq) 



T 



[1 



r' C r'-r 



<50(r) <50(r') 



dr' 



(2.8) 



follows which is analogue of equation (2.14) in [9] and determines equilibrium DF. 
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Unfortunately, to the best of my knowledge, non- stationary solutions to "generating 
equations" like (12.51) (or (7.9) from [9]) never were investigated by Bogolyubov or any 
other authors. However, past experience in the BBGKY equations (see e.g. pi)]) points 
to desirability of a change of variables, i.e. transition from DF to some suitably defined 
"correlation functions" . With this purpose, let us discuss hypothetical equalities 

F n (t, R, rW, P, p(»)| R ; u) = F (t, R, P| R ; u) (r W I R; v) ]$ =1 G m ( Pj ) 

or, equivalently, F{t, R, P, tp | R ; v\ = F (t, R, P| R ; u) ^ eq) {<f) | R; u} (recollect 
that 0(r) = J ip(r, p) G m (p) dp). They state that correlations of atoms with wan- 
dering BP always stay the same as with pinned BP. It seems reasonable in view of 
thermodynamically equilibrium character of the wandering. Nevertheless, it may be 
true only if all possible BP's positions are statistically equivalent. The latter in our 
case is not true since translation symmetry is destroyed by information about BP's 
start point Ro . Correspondingly, above equalities are incompatible with equations 
(12. 2p . For example, substitution of equality for F\ to "collision integral" in equation 
for Fq turns it into zero as if BP does not interact with atoms at all. 

The aforesaid shows that, first, BP's wandering produces specific non-equilibrium (in 
statistical sense) "historical" correlations between its total path AR(t) = R(i) — Ro 
and current state of surrounding medium. Second, we can adequately separate these 
correlations from equilibrium ones if define them as follows: 

F{t, R, P, i) | R ; v} = V{t, R, P, tj; | R ; v} F {eq) {(j) \ R; v} , (2.9) 
V{t,R,P,^|R ;z/} = Vq /- / V " II^ r i'Pi) ' 

n=l U - Jl Jn j=l 

where V n = V n (t, R, rW, P, p (n) | R ; v) are corresponding correlation functions (CF) 
and V their GF. In particular, evidently, Vo(t, R, P| R ; v) = F (t, R, P| Rq; v) and 
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F 1 (t,R,T 1 ,P,p 1 \R ;u) = 

= F (t,R,P|R ;z/)F 1 (e9) (r 1 |R;z/)G m (p 1 ) + Vi(t, R, r 1; P, p x | R ; v) 1 ' > 

In terms of CF the initial conditions (12. 3p and "normalization to volume" conditions 
("weakening of correlations" at infinity) take very simple form: 

V n (t = 0) = 5„, 5(R-Ro) , V{0, R,P, V|Ro;^} = 5(R-Ro) , (2.11) 

V n>0 (t,... r k -> oo ...) -> (2.12) 
Substituting (12.91) into (12.51) one obtains equation for GF of "historical" correlations: 

t'fa'jj) = {/[l+0(r)]^(R-r)^{r,0|R^}dr}^ + ^ + 

+ ^ [1 + V(zi)][l + V(za)] $ a (ri-r 2 ), z/C{r 2 ,0|R; z/} G m (p 2 ) — ^-y 

where new functional C is defined as 

„, , .„ . tflnj^MlR; z/j 

C{r, R; u} = IJ\ ] = 2.14 

= Fi eq \r\R;v) + J2— - C n+1 (r, r x ... r n | R; v) J[ 0(r,) drj 

n = l n - J J j=\ 

Complication of equation (12.131) in comparison with (12.51) is pay for simple conditions 
f l2.lip - fj2.12p . Corresponding equations for CF also are more complicated than ( 12. 2ft . 
Therefore here we write out them only for extreme but interesting case of "BP in ideal 
gas" (when $ tt (r) = , i.e. atoms do not interact with themselves): 

dV P dV d 
+ v 



dt M dR dP 



9V n>0 d 
— - — = \H n , V n + v — — 
dt [ ni nJ OP 



y"$' 6 (R- ri )T4 , 

/ & b (R-r n+1 )V n+1 + (2.15) 

Jn+l 



n f P 8 \ 

+ TJ2 Hj,n) G m (Pn) E'(v n - R) ( — + — J V n . 



i j 
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where H n = P 2 /2M +£"=i [ p 2 /2m+$ b (R - r,) ] , E{r) = exp [ - $(r)/T] , £'(r) = 
V.B'(r) = — <&'(r)E(r)/T and V(j,n) denotes operation of transposition of arguments 
2Tj = {rj,pj} and x n = {r n ,p n }. Thus natively bidiagonal BBGKY hierarchy in 
terms of CF becomes tridiagonal. 



3. Invariance group of equilibrium generating functional 

Equilibrium CF C n in (I2.14p are determined by first equality in (12.141) : 

C a (r,n|R;i/) = F^\v,v l \^v)-Fi eq \v\^v)F[ eq \v x \R ] v) , 
C 8 (r,n,r 2 ) = F 3 (e9) (r, ri ,r 2 ) + 2^)^^)^)^)- 

_ Fl ^)( r)F N)( ri;r2 ) _ F^Vi)if 9 V,r 2 ) - F^^F^Cr.n) , 

and so on. In the second of these expressions for brevity we omitted arguments R 

and v . The conditions of normalization to volume, i.e. weakening of correlations at 

infinity, which in essence establish existence of the thermodynamical limit [9j [11] , do 

mean that all these CF vanish when distance between any two atoms goes to infinity 

and turn into usual equilibrium CF when BP is moved off to infinity: 

(7 n+1 (r,r 1 ...rjR;i/) -»• at r 5 - r -> oo , , g ^ 

C n (ri ...r n | R;z/) -> C n {r x ...r n ; v) at R - r, -> oo 

Moreover, under sufficiently short-range interactions all these limits are achieved in a 

fast (absolutely) integrable way, which will be assumed below. 

Next, consider equation (12.81) rewriting it in the form 



dr T 



C{r, <t> |R; i/} = - J[ 1 + 0(r') ] - r) " } dr' + 

+ C{r, |R; u} £ J [ 1 + 0(r') ] <^(r' - r) C{r', |R; „} dr' (3.2) 

as equation for the functional C . The latter interests us here only so far as it influences 
equation (12.131) for GF of "historical correlations". In this respect, we have (perhaps, 
for the first time) to pay attention to some important properties of this GF. Firstly, 
due to integrability of the asymptotic (13.11) . C can be extended to bounded functions 



<p(r) which do not turn to zero at infinity, particularly, to constants, and introduce 
objects as follow: 



C(a,v)= lim C{r,a|R;z/} = 1 + V 



n = l 



TV. 



1 Jn 



C n+ i(r,ri ...r n ; v) , (3.3) 



C CT {r,0|R; v) 



C{r, a + (f> |R; v} 



(3.4) 



where a = const , J ... = J ... dr n , and integrals in (13.31) are factually independent 
on r . Secondly, replacement 0(r) — > o + 0(r) in (13. 2p and elementary algebraic 
manipulations transform (I3.2p to equation for C a : 



d $1(t - R) 



dr 



T 



+ C a {r,(p\R;u} 



C a {r,(j)\R; v} = - 



uC(a, u){l + a) 
T 



1 + 

(h 

1 + 



l + a 
r 



l + a 



al j 5[0(r)/(l + ( r)] 
^(r'-^CJr'^lR;//}^' (3.5) 



It differs from (13.21) only by scale transformation of the functional argument, 0(r) 
0(r)/(l + a) , and replacement of the density v by 



f(cr, v) = uC(a, v) (1 + a) 



(3.6) 



Third, formal solution to (13.21) in the form of series (12. 14j) is unambiguously determined 
by the normalization conditions, that is conditions (13.11) plus quite obvious equality 
lim R_ r ^oo C{r, <p = |R; u} = 1 which determines first term of the series. Fourth, 
functional C a defined by (I3.3p - (I3.4I) satisfies same conditions since (13. ip imply 



— / ... / C n+fc+ i(r,ri...r n+fc |R;i/) -> at r,- 



(3.7) 



at least, if one understands (13.11) in the sense of absolute integrability of CF C n+ i and 
speaks about not too large values of the density t . 

Summarizing all that, we can conclude that solution to equation (13.51) is nothing 
but C a {r, (ft |R, u} = C{r,0/(1 + a) \ R; v(a, u)} . With taking into account (13.41) 



2 An useful information in this respect comes from rigorous consideration of "group properties of 
correlation functions" ITT]. 
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and (I3.6P this means that at arbitrary (admissible in definite sense) constant a and 
bounded function (p = 0(r) equality 

vC{r,a + <f>\R;v} = C (r, -A_ | R; v(a, u) ) (3.8) 

1 + cr [ 1 + a J 

is valid . It can be rewritten also as 

f(a) C{r, | R; v} = C{a, v) C jr, - 1 1 R; v(a, v) j = C{r, | R; v) , (3.9) 

where functions C(a, v) and ^(cr, v) are connected by means of ( 13. 61) . and the left 
equality defines one-parameter family of such transformations of arguments of func- 
tional C{r, cf) | R; 7v} which do not change its value. It is not hard to verify that this 
is group described by the composition rules 

f(a 2 )f(a 1 ) = f(a 1 + a 2 + (J 1 a 2 ) , 
v(a 2 , v(ai , v)) = v(ai + a 2 + a x a 2 , v) , (3.10) 
C{p 2 , (1 + (Ti) C{o\ , v) v) C(cri , v) = C(<Ji + a 2 + o-icr 2 , ^) 

with restrictions a > — 1 , 0(r) > — 1 . The latter are clear in the light of that 

t;{r| 0, R; i/} = [1 + 0(r)] ^' ' 1 = z/ [1 + 0(r)] C{r, | R; u} 

represents mean concentration of atoms under external potential U(r) which is related 
to 0(r) by 0(r) = exp [ — U{r)/T] — 1 (see e.g. [12J). Substitution a = exp (a) — 1 
gives T(a 2 ) T(a\) = T(a\ + a 2 ) thus eliminating the restrictions. 

Infinitesimal form of (13. 8ft or (13.91) looks best if written through particular CF : 

| x ( z/ ) + [l + x ( z/ )] z/ |_| Fi ^)( r | R;z/ ) = ujc 2 (r,v'\R;u)dr' , 
jnx(z/) + [ 1 + k{v) } v C n (ri ... r n | R; v) = v J C n+1 (r 1 ... r n , r'|R; v) dr' , 



X 77 



dC(a, v) 
da 



J cr=0 



v J C 2 (r,0; is)dr (3.11) 
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The function x(u) is known (see e.g. p3J) to be directly related to a state equation 
of the system: 1 + x(v) = T(di>/dV)T , where V denotes the pressure. Notice that 
in the framework of the grand canonical ensemble substantially similar relations can 
be easy derived by differentiation of DF in respect to the activity. 

4. Invariance group of generating functional 
of historical correlations 

Now, let us show that solution to equation (12.131) has invariance properties quite 
similar to (13. 9p . Since initial condition to this equation (see ( 12. lip ) does not depend 
on variables ip = ip(r, p) and v at all, solution to it is completely determined by 
structure of operators £ and £ ' and conditions (12.121) . The latter formally allow to 
extend functional V{t, R, P, ip |Ro; u} (like C before) to arguments a + ip(r, p) , with 
a = const , in place of ^(r, p) . The fact that limit in (12.12ft is achieved fast enough 
to indeed ensure this extension can be confirmed afterwards. Besides, thanks to (12.121) 
variable ip(x\) in expression £V inside ( 12.131) (see definition ( 12.61) of the operator 

£ ) can be shifted by arbitrary constant: 

f pi 8 5V f Pl d 5V 

J 1 m or i Qip{xx) J x m Ori bip{x\) 

where a = const , for instance, a = 1 . This is important difference of £V from 
£ T . Consequently, taking in mind action of £ onto V , one can write 

Further, let us carefully consider operator £' (see also (I2.13P ). In contrast to £ , it 
depends on the density v . Nevertheless, with the help of equality (13. 8p it is easy to 
make sure that it obeys the same relation if transformation of argument ip(r, p) is 
accompanied by transformation of argument v in accordance with (13. 6p and f 1 3 . 3 [) : 

C'^t^)^' 1 ^,^, ,,^,!! ) (4-2) 
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Formulas (14 .11) and H4.2[) just imply the noted invariance property of solutions of fl 2 . 1 3 p : 

V{t, R, P, a + i> | R ; v} = V jt, R, P, | R ; v(a, v) J (4.3) 

or, equivalently and similarly to (13. 9p . 

f(a) V{t, R,P, ^ | R ; u} = V{t, R, P, ~ 1 1 R o! u(<7, i/)} = 

1 + a 

= V{t,R,P, V|R ;z/} , (4.4) 

where left inequality together with ( 13.61) and ( 13. 3D defines action of the above described 
group onto GF of historical correlations. Expansion of H4.3[) into series over ip yields 

V (t, R, P| R ; v(a, v)) = V {t, R, P| R ; v) + (4.5) 

n = l n! ^ 

l>(<7, I/) 



2/ 



1 + <T)*> 



V fc (t,R,r( fc ),P,p( fc )|R ; u(<t,i/)) = Vfe(i, R, r (fc) , P, p (fc) |R ; i/) + 

+ E^-/ -/ ^ +n (t,R,r( fc+ "),P,p( fc+ ")|R ;z/) (4.6) 



n! 

n = 1 



'fc+1 ./fc+n 

Corresponding infinitesimal (in respect to a ) relations are similar to (13.111) : 

|nx(i/) + [l + x(z/)]z/|-| K(t,R,r (n) ,P,p (n) |R ; i/) = (4.7) 

= i/ / K+i(t,R,r (n+1) ,P,p (n+1) |R ; v) 

Formulas (14 . 5 [) - ( 1476T) can be interpreted as "virial expansions" of probabilistic law 
of BP's random wandering and historical correlations between BP and medium, with 
those difference from usual virial expansions of thermodynamic quantities [13] or kinetic 
coefficients [H] that here decrements of the density do figure instead of its full value. 
However, in the limit v — > , cr^oo, vo = const these relations take quite usual 
form. The simplest of them, (14.51) . recently was obtained [121 EE EH] in other way, 
starting from the "generalized fluctuation-dissipation relations" [XTf [18]. 
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As far as I know, exact relations of this kind never before were under consideration. 
It would be rather hard to extract equalities (14.71) (all the more, (I4.5p - (I4.6I) ) directly 
from the BBGKY equations (12. 2\\ . This is possible only in the case of "BP in ideal gas" 
(the reader can try to see how (14.71) follows directly from equations (12.151) for CF). 

5. On principal consequences from the virial relations 

All the exact "virial relations" (I4.5I) -( 14T71) are automatically satisfied by exact solution 
to BBGKY equations. Therefore these relations can be applied to "testing of statistical 
hypotheses" about the solution or constructing approximations to it. Moreover, in the 
just mentioned limit v — > , er^oo, va = v' = const, with the help of formulas 
(13. 3p . (13.41) and (13.61) . it is easy to transform equalities (I4.5l) - (l4.6p into explicit formal 
representation of exact solution of BBGKY hierarchy for V n (t, R, r", P, p( n )|Ro; v') 
as a power series in respect to v' . 

It should be emphasized that (I4.5p - (I4.6I) connect random walks of BP in two media 
whose densities may form arbitrary large or arbitrary small ratio, v(a, v)jv . Therefore, 
any cutting off the series in (I4.5p - (I4.6I) would present not an approximate but incorrect 
result. This means that from the viewpoint of exact theory all historical correlations 
always are equally significant. Even in the "Boltzmann-Grad limit" ( "dilute gas limit" ) 
when gas parameters tend to zero, 47rr^z//3 — > 0, 47rr^z//3 — > (with r a ^ standing 
for radii of interactions of atoms and BP) while free paths A;, = (7rr^z/) _1 and A a = 
(7rr^z/) -1 stay fixed, or in the limit of ideal gas (where A a = oo but A& is finite). In 
both cases, formulas (I4.5p - (I4.7I) simplify to 

■../ W,R, r <«+W+«>|R^) (5.1) 

Jn+l Jn+k 

We must conclude that Boltzmann's kinetics which does not know correlations (except 
may be two-particle one) is not true "zero-order approximation" in respect do density. 
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As an illustration, let us consider BP in ideal gas and apply relations (15.11) to test 
"statistical hypothesis" that a correct approximate solution of exact equations (12.151) 
can be obtained under neglect of three-particle and higher correlations. Such hypoth- 
esis always is (directly or indirectly) involved into derivation of kinetic equations for 
"sufficiently rarefied" gas. Examples can be found e.g. in [91 [TOj [HI HH [20] Q Set- 
ting V2 = in the second ( n — 1 ) of equations (12.151) and then inserting result of 
its integration to the first one ( n = ) we come to closed Boltzmann-Lorentz kinetic 
equation [19], [20] for Vo(£, R, P|Ro; v) ■ Asymptotic of solution to this equation at 
t ^> t = Ab/i>o (with i>o ~ a/T/M ) is wittingly Gaussian: 

V (t, AR; u)=J V (t, R, P|R ; u) dP -> ^['^(^ > ( 5 - 2 ) 

where D is BP's diffusivity, D = v A b oc u^ 1 . Undoubtedly, this is a complicated 
function of the density v . At the same time, in the view of exact relations (15.11) the 
statement V2 = implies that Vo(t, AR; v) should be purely linear function of v ! 

So strong discrepancy prompts that our hypothesis is erroneous, asymptotic (I5.2p is 
doubtful, and thus we should return to BBGKY equations. More correct approach to 
their approximate solving, with including correlations of any order, was suggested in 
[8] (or see [22] ) and developed in [23] . It confirmed the guess [24] that actual molecular 
random walk represents a diffusive random process possessing scaleless "flicker" (i.e. 
with 1/f -type spectrum [7]) fluctuations of diffusivity (as well as mobility) of BP. 
In corresponding asymptotic of Vo(t, AR; v) , in contrary to (I5.2p . the exponential is 
replaced by a function with power-law long tails (cut off at distances ~ v§t ) [23] : 



V (t, AR; v) 



r(7/2) 



AR 2 



-7/2 



[AnDt] 3 / 2 [ ADt 
where 0(x) ~ 1 at x < 1 and 9(x) — > at x > 1 in a very fast way. 



e|£S] , (5.3) 



3 The work [5D] on gas of hard spheres gave example of illusory discard of correlations, as it is 
discussed in [5T] . 
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Origin of the diffusivity fluctuations is trivial. Heuristically, that is indifference of 
the system to a number and relative frequency of BP's collisions with atoms (all the 
more, to distribution of collisions over impact parameter values) [HI E2]. At that, 
higher historical correlations described by V n> \ are caused by complicity of particles 
in (uncontrolled and therefore scaleless and may be non-ergodic) fluctuations of relative 
frequency of BP's collisions (see also [121 ESI ED E3])- 

Let us demonstrate that one can come to analogous conclusions without solving 
BBGKY equations but instead resting upon their invariance group and besides the 
trivial fact that all the DF F n always are non-negative. 

From non-negativity of F\ and identity (12.101) we have 

V (t, AR; v) [ Fi eq \r\R;u)dr + f ^R^Ro; v) dr > , (5.4) 
Jq Jq 

where V\(t, R, r|R ; v) = J J Vi(t, R, r, P, p|R ; v) rfpdP , and Q is any region in 
the space of vectors r — R . Introduce 0(5, t, AR; v) = Q(5) to be minimum (in the 
sense of volume) of all regions Q what satisfy condition 

(5.5) 



1 V x dr - [v 1 dr 


< $ 


J V x dr 


Jq J 







where < 5 < 1 . One can easy justify that (15.41) and ( 15. 51) together imply inequality 

V {t, AR; v) I F[ eq \r\^v)dx + {1-5) [ T4(t,R,r|Ro; u)dr > (5.6) 

Combining it with exact virial relation (14 .7p for n = and taking into account the 
equality 1 + x{y) = T [dv /dV)r (see Section 3) after clear reasonings one obtains 

T V) + m, AR; v) V Q (t, AR; u) > , (5.7) 

where quantity Q(t, AR; v) is defined by 

m,AK;u) = min — ^ I F^(r| R; u) dr w min ^ll^2 v ) ( 5 . 8 ) 

0<5<1 1 — 6 J Mg) 0<8<1 1 — 

According to the aforesaid, if Vi(t, R, r|R ; v) as a function of r (at given other ar- 
guments) has constant sign, then Q(t, AR; v) represents volume occupied by the pair 
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correlation or, briefly, "pair correlation volume" . Otherwise £l(t, AR; v) is something 
smaller than this volume. 

Now discuss hypothetical asymptotic (15.21) from the viewpoint of inequality (|5.7p . It 
says that the hypothesis can be true only if 

7^, x dln/J /AR 2 3\ AR 2 3 . , 

uQ(t, AR; v) > -vT—— — J -> — , 5.9 

v ' ' ; ~ &P \ ADt 2 / ADt 2 ' v ' 

where last expression concerns dilute gas. In other words, if quantity vQ(t, AR; z/) , 
all the more the pair correlation volume, measured in units of specific volume per one 
atom, l/u , is not bounded above. In opposite, if it is bounded above, 

i/fi(t,AR; v) < c x = const , (5.10) 

then Gaussian asymptotic (15.21) is forbidden. Instead, inequality (15. 7p allows for a 
generalized diffusion law as, for instance, 

V (t, AR; v) - / 1 — * ( 6 ( , (5.11) 



^— j^- + a^f(z) > , a = - + a 



vT 



(5.12) 



where function *&(z) should satisfy inequality 

d\nD~ 
dV 

(supposing that D falls when pressure grows). Consequently, ^{z — > oo) oc l/z a . 
In the case of gas, a = 3/2 + c\ . Formula (15. 3p corresponds to ci = vQ, = 2 . 

Hence, the theory inevitably leads to statistical correlations which have unlimited 
extension either in space, as under variant ( 15. 9ft . or in time, as under alternative variant 
(15.101) . It remains to ascertain what of the variants is closer to exact solution of BBGKY 
equations. From physical point of view, the second one certainly is preferable, at least 
if speak about gas. Indeed, as it is clear from equations ( 12. 15ft . a source of correlations 
between BP and atoms is their collisions. A collision realizes at such disposition of 
BP and atom when vector r — R lies in the "collision cylinder" which is oriented in 
parallel to the relative velocity p/m — P/M and has radius ~ . At that distance 



17 

between colliding particles should not be much greater than A = min (A OJ A&) , since 
otherwise their collision almost surely will be prevented by an encounter with the rest 
of atoms. Consequently, at any values of momenta (all the more, at any t and AR ) 
it is natural to estimate the pair correlation volume as Q ~ 2A -kt\ . If BP is mere 
marked atom |23j . then A = A& and Q ~ 2Af,7rr^ = 2/u . As the result, we arrive 
to CT)} with ci » 2 . 

6. Conclusion 

Thus, to resume, the problem about thermodynamically equilibrium random walk of 
a molecular "Brownian" particle (BP) in a fluid was formulated as a problem of classical 
statistical mechanics of a system of (infinitely) many particles. Corresponding BBGKY 
equations for distribution functions (DF) and equation for their generating functional 
were considered in terms of correlation functions (CF) introduced so that they extract 
statistical correlations between total path of BP during all its observation time and 
current state of the medium ("historical correlations"). We showed (in case of Gibbs 
canonical ensemble) that generating functionals of both equilibrium DF and time- 
dependent CF are invariant in respect to definite continuous group of transformations 
of their arguments including density of the medium (mean concentration of particles). 

The found invariance group produces a sequence of exact "virial relations" which 
connect full sets of DF or CF taken at different values of the density and therefore can 
serve as quality test of approximate solutions of the BBGKY equations. In this respect, 
we demonstrated that conventional "Boltzmannian" approximation to kinetics of BP 
in dilute gas (or ideal gas) is incorrect. The matter is that it rejects correlations of third 
and higher orders, while actually correlations of any order are equally important, even 
(and most of all) under the Boltzmann-Grad limit. The virial relations automatically 
allow for all the correlations and therefore imply significant restrictions on possible 
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forms of probabilistic distribution of the BP's path. In particular, they rather surely 
(especially in case of dilute gas) forbid Gaussian asymptotic of this distribution. What 
is possible instead of it is an automodel diffusive asymptotic possessing power-law 
long tails (cut off at ballistic flight length). This conclusion well agrees with previous 
approximate solutions of BBGKY equations [H [221 123]. It means that in real gas 
random BP's trajectories are so much unique that, in contrast to the Lorentz gas, they 
can not be divided into statistically independent constituent parts. In other words, 
they can not be imitated by "dice tosses" . In this sense actual molecular random walk 
is non-ergodic M, EH [24] . 

In should be noted that results of this work can be easily extended to thermody- 
namically non-equilibrium walk under influence of an external force applied to BP from 
start of its observation (in this respect see [T2l [16]. 122]). From the other hand, it would 
be interesting to generalize the mentioned invariance group to more usual problems 
where not some select particles but hydrodynamic fields are in the centre of attention. 

I would like to acknowledge my colleagues from DonPTI NANU Dr. I. Krasnyuk and 
Dr. Yu. Medevedev for many useful discussions. 
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